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Abstract 

We study gravitational waves in viable f(R) theories under a non-zero background curvature. 
In general, an f(R) theory contains an extra scalar degree of freedom corresponding to a massive 
scalar mode of gravitational wave. For viable f(R) models, since there always exits a de-Sitter point 
where the background curvature in vacuum is non-zero, the mass squared of the scalar mode of 
gravitational wave is about the de-Sitter point curvature ~ W~ 6e eV 2 . We illustrate our results 
in two types of viable f{R) models: the exponential gravity and Starobinsky models. In both cases, 
the mass will be in the order of 10 _33 ey when it propagates in vacuum. However, in the presence of 
matter density in galaxy, the scalar mode can be heavy. Explicitly, in the exponential gravity model, 
the mass becomes almost infinity, implying the disappearance of the scalar mode of gravitational 
wave, while the Starobinsky model gives the lowest mass around 10 -24 eV, corresponding to the 
lowest frequency of 10~ 9 Hz, which may be detected by the current and future gravitational wave 
probes, such as LISA and ASLROD-GW. 

PACS numbers: 04.30.-w,04.50.Kd 



* louis.lineage@msa.hinet.net 
f |g9522545@oz.nthu.edi!tw| 

* geng@p hys.nthu.edu. tw| 



1 



I. INTRODUCTION 



The accelerating expansion of our universe has been supported 
observations, such as typeja supernovae jl, 2|, large scale structure 



background radiation |5|, |6( and weak lensing Q . In order to explain this late time acceler- 
ation one can either introduce dark energy, a new form of matter, or modify Einstein's 
general relativity, i.e., the modification of gravity. One simple way to modify general rela- 
tivity is to promote the Ricci scalar R in the Einstein-Hilbert action into an f(R) function, 
which is the so-called f(R) theory 9Ml3|. A viable model of f(R) can generate a late-time 
accelerating expansion of our universe, have the radiation-dominated stage followed by the 



jy various cosmological 
|, |4| , cosmic microwave 



15|, and be consistent with the solar-system constraint under 
23|. 



matter-dominated one Jl 
chameleon mechanism [16 

In 24|, Chiba showed that an f(R) model will allow a new scalar degree of freedom. This 



corresponds to a new scalar mode of gravitational wave besides the ordinary tensor one of 
general relativity. This new scalar mode will be massive and propagate as a longitudinal po- 
larization. Various discussions and predictions about this extra scalar mode of gravitational 
wave have been given in the literature 25N32|. However, most of them were concentrated 



on either quadratic or inverse-curvature type of f{R) models, which is highly restricted by 



12 



20 



33 



36|. In the article, we will study gravitational wave in 



the observational results 
viable f(R) models. 

The mass of the scalar mode of gravitational wave in viable f{R) models could be quite 
different from quadratic and inverse-curvature ones. In vacuum, it will be of the order of 
the Hubble constant because all viable f(R) models need to have de Sitter points which 
have a non-zero background curvature about square of the Hubble constant. However, 
when gravitational wave propagates in the galaxy region, the local density of dark matter 
and baryonic matter will contribute a tremendous background curvature although it is still 
much smaller then curvature generated by star. This background curvature might make 
some viable f{R) models return to the ordinary GR very fast. Therefore, the scalar mode 
will become extremely massive in these cases and, hence, prevent the observable propagation 
of the fifth force in our galaxy and solar system. On the other hand, viable f(R) theories 
pass the solar-system constraint by means of chameleon mechanism, which do not require 



very heavy scalar modes by using the thin-shell argument 



2l| . In these models, the 



scalar modes may still be detectable. 

This paper is organized as follows. In Sec. [TXJ, we review the field equations and lin- 
earizations in f(R) theories and demonstrate the modifications of gravitational waves. We 
apply the analysis on two explicit viable f(R) models: exponential gravity and Starobinsky 
ones in Sec. IIHI The results and discussions on the scalar modes of gravitational waves in 
the inner galaxy region is presented in Sec. IIVI Conclusions are given in Sec. |VJ We use 
units of ks — c — H — 1 and the gravitational constant G = Mpf with the Planck mass of 
Mpi = 1.22089 x 10 19 GeV. 

II. GRAVITATIONAL WAVES IN f(R) THEORY 
A. f(R) Gravity 

We start by considering a general Einstein-Hilbert action 

S = ^J d A x^-gf\R) + S m {g^, (2-1) 

where f(R) is an arbitrary function of the Ricci scalar R, S m is the action of the matter 
part and k 2 = 8irG. In the metric formalism, we vary the action (12. ip with respect to g^, 
and the modified Einstein field equation can be obtained as 

f'(R)Rnu - \f(R)g»» + (g^n - v„V„) f'(R) = k% u , (2.2) 

where a prime denotes the derivative with respect to R, V M is the covariant derivative and 
□ = g^V^Vu is the d'Alembert operator. The trace of the field equation ( 12. 2ft gives 

f'(R)R - 2f(R) + 3D f'(R) = k 2 T, (2.3) 

where T = g^T^ = —p + 3a 2 P is the trace of the matter energy-momentum tensor, and a 
is the scale factor. 

For f(R), the de Sitter stage is a vacuum solution with a positive constant background 
curvature Rd, which is assumed to be homogeneous and static. Consequently, one has 

V ll f(R d )=0 and f'(R d )R d = 2f(R d ). (2.4) 
Moreover, from Eq. (12. 2p . the Ricci tensor satisfies -R MV |_R d = g^vRd/^- 
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B. The Weak-field Approximation 



In order to investigate gravitational wave in f(R) theories, we need to study the linearized 
theory of f(R) gravity. Consider a small perturbation from the FRW metric: 

sv = v + V" ( 2 - 5 ) 

where \hfi V \ <§C 1 is the perturbation and ~g^ u = diag(—l,a 2 ,a 2 ,a 2 ) is the FRW background 
metric. If the evolution of the system is much shorter than Hubble time, we can approximate 
the background spacetime to be nearly the Minkowski one with g « rj^ = diag(—l, 1,1,1). 
We keep the theory to be the first order in and neglect terms higher than O {h 2 ). Thus, 
the inverse of the metric tensor is given by 

gt* = gt™ _ (2.6) 

Note that all indices are raising and lowering by the background metric g . In the metric 
formalism, the perturbation of connection is 

SFLp = (dph ap + d a h^ - dfjicp - 2h ltv T a p) , (2.7) 
where T a a is the unperturbed connection. The only non-vanishing components of T a/3 are 

f *„ = r 0j = H5) and r°- = a 2 H8 ij , (2.8) 

where H = a /a is the Hubble constant. The deviation of the Ricci tensor from the back- 
ground curvature is j^l 

5R afj = d,5T» - dpsr* + r^sr» + r^sv^ - r»sr» - f* ar^- (2.9) 



C. The Scalar Mode h 



f 



The different between gravitational waves in f(R) and general relativity is that it contains 
an extra scalar degree of freedom in f(R). This comes from the non- vanishing trace of the 
field equation. Eg . (|2.3I) can be viewed as equation of motion for a scalar field $. By the 
identifications 



38 



39 



K ' d$ 3 



(2.10) 



we obtain the Klein-Gordon equation for the scalar field $: 

dVeff 



(2.11) 



In order to have a stable perturbation of spacetime, we must require the background scalar 
$o t° stay at the stable minimum of the effective potential V e //; i-e., 

dV ef f 



and 



d$ 2 







> 0. 



(2.12) 



(2.13) 



In vacuum, Eq. (12 . 1 2[) just gives us the condition for the de-Sitter point curvature (12.41) . 
while Eq. (12 . 1 3[) requires the mass of the scalar mode to be positive. 

Perturbing the trace of the field equation ( 12. 31) with a nonzero constant background 
curvature Rq which satisfies Eq. (12 . 12[) yields 



3Q5f + R 5f + f'(Ro)5R - 25 f = 0. 



(2-14) 



Using the relations S 
for the scalar mode 



27 



Rq)5R and 5f = f"(Ro)5R, we obtain the massive wave equation 

(2.15) 



mjx 



where hf = 8f'/f'(Ro) is the field of the scalar mode and 



m 



f'(Ro) 



Ri 



(2.16) 



3 \f"(Ro) 

is the mass squared of it. Note that m? s = V"ff(<&) 19|. For any viable f(R) model, the 



condition m 2 s > is needed for the stabilit y of the cosmological perturbation and to prevent 
the field from being a tachyon 
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40 



4JJ. 



For the FRW metric, Eq. (12. 15ft should be expressed as 

-81 + ^ - 3Hd ) h f = m 2 s h 



(2.17) 



where the term —3Hdo gives a damping factor caused by the expansion of the universe. To 
illustrate the solution of Eq. (I2.17p . we take the de Sitter universe with a constant H. In 
this case, the solution is a damped plane wave 



A(k)e 2 exp (iq^Xa) 



(2.18) 



where g M = (cu m ,k), u m = \Jk 2 /a 2 + m 2 — \H 2 is the angular frequency and A(k) is the 
amplitude. For simplicity and without loss of generality, we take a = 1 and neglect the 
damping effect as k 2 /a 2 3> H 2 . As a result, Eq. ( 12. 15ft leads to a simple plane wave solution 



A(p)exp(iq ,1 x ll ), (2.19) 



with uj m = y k 2 + m 2 . We can see that m s is the cutoff frequency of the scalar mode of 
gravitational wave. For u m < m s , the wave vector becomes imaginary. The waveform is an 
exponential decay in distance, i.e., hf oc exp(— k ■ x). Thus, the scalar will not propagate in 
space below the cutoff frequency. The massive scalar mode will not propagate at the speed 
of light with the group-velocity 



(2.20) 



D. The Tensor Mode hT„ 



I w 



Perturbing the field equation (12.21) under the de-Sitter curvature R d leads to 

f'(R d )5R, u + R^UJf - l --g^Jf - \K v f(R d ) + (g^D - V M V„) Sf = 0. (2.21) 
Using 5f = pf, R^\ 

R d = 9uuRd/^ and the condition for the de-Sitter stage curvature 
(glD , Eq. (I2T2T]) becomes 

SR^ + \g^R d h f - \gjjhj - + (V D - V mV.) h f = 0. (2.22) 

The total perturbation of metric can be decomposed into the tensor part h T and scalar part 



If 



in the way that 

V = h lv + h%> ( 2 - 23 ) 
where h^ u = bg^hf and b is an unknown factor which will be determined later. Note 
that we require the tensor mode fvF v to be traceless because we do not want it to give any 
contribution to the purturbation of Ricci scalar 5R and couple to the scalar mode hf. Since 
the theory is linearized, the perturbation of Ricci curvature tensor can also be seperated 
into two parts 

6R IJV = SR% / + SR^, (2.24) 
(3 



where 5R? U and SR^ V represent the perturbations contributed from the tensor mode h^ u and 



scalar mode h^ u , respectively. Jn the de Sitter universe, we have 8qH = and R, 
Thus, SR^ U can be written as 



"fill 

12H 2 . 



*R%, = ~b 



Inserting this into Eq. (I2.22p . we obtain 



g^U + V^V, ) h f . 



1 



1 - - ) g^D - (1 + b) V M V, + - (1 - 6) g^R 



I- r_ 



(2.25) 



= 0. 
(2.26) 



Since the wave equation for the scalar mode ( I2.15P does not involve any off-diagonal term, 
and we do not want the coupling between scalar and tensor modes, the term V^V^/i/ should 
not appear here. To cancel the V^Vj, term, we must pick b = — 1 
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5R 



i"' 



9 



3D + R d 



L 
f" 



h f = 0, 



3lfl . This gives 

(2.27) 



where the terms in parentheses vanish by using the wave equation for the scalar mode 
f)2.15p . It is clear that the extra scalar degree of freedom in f(R) can totally decouple from 
the ordinary tensor mode: 



6Rf 



0. 



"nv (2.28) 
Similar to gravitational wave in GR, we assume the tensor mode to be divergenceless, 



which means that it has to satisfy the Lorenz gauge condition, i.e., d^h^ = 



27j, 



3l|, and 



be transverse, hj^ 



h,Q = 0. The perturbation of Ricci tensor for the tensor mode can be 



simplified in this 



SR 1 = 

p,v 2 



-d 2 + % + Hd 

or 



AH 2 h 



(2.29) 



We further define hjj = a 2 l-iij to absorb the effect of expansion of the universe 42], and 
choose the tensor mode to propagate in the z direction with 





(k 


h x 


°1 








-h+ 





(2.30) 




1° 





°) 





where h + and h x are the plus and cross polarizations, respectively. The first and second 
time derivatives of hf^ can be expressed as 



a 2 (d + 2H) Uij and d 2 hjj = a 2 (d 2 + 4Hd + AH 2 ) U 
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i.r 



(2.31) 



Inserting these and Eq. (I2.29P into (I2.28p . we obtain the wave equation for the tensor mode 

(2.32) 



a 2 



-d 2 + ^ - 3Hd ) Hij = 



or 



UK 



0. 



a 



+, x. 



(2.33) 



Note that Eq. fl2T32l) is equivalent to Eq. (5.1.53) in Ref. |37( and Eq. (5.61) in Ref. (42|. 
Clearly, the tensor mode is exactly the same as that in GR when a traceless gravitational 
wave propagates in a non-zero de-Sitter curvature Rd background. In what follows, we will 
concentrate on the scalar mode of gravitational wave. 



III. GRAVITATIONAL WAVES IN VIABLE f(R) MODELS 



The conditions for a cosmological viable f(R) model include (i) the positivityof the ef- 



fective gravitational coupling, (ii) the stability of 



cosmological perturbations 



(iii) the stability of the late-time de-Sitter point 14|, 44 ^46 [ , (iv) the asymptotic behavior to 



36 



40 



41 



ACDM at the high curvature regime, (v) the solar system constraint, and (vi) the constraint 



from the violation of the equivalence principle |18l . |19|. \21\ . The typical examples of the 
viable f(R) models are Hu-Sawicki 19 1, Starobinsky 36 1, Tsujikawa J22I and exponential 
gravity models {47 1. To illustrate our results, we will concentrate on the exponential and 
Starobinsky models. Our study can be easily extended to other viable models. 

Usually a viable f(R) model can use chameleon mechanism to pass the constraints from 



2l|. In this light-mass scalar is 



the solar-system and equivalence principle 
allowed by introducing the thin shell condition. As a result, the mass of the scalar mode 
does not have to be very heavy 18] . 



A. ACDM 

We can take the cosmological constant model or the ACDM model as a special case of 
f(R) with 

f(R) = R-2A, (3.1) 



where A is the cosmological constant. 



In this case, the mass of the scalar mode is infinite, i.e., m? s = oo, which requires infinite 
large energy to excite the scalar mode. Clearly, there is no scalar mode in the ACDM model. 
Although the de Sitter curvature R d is not zero, i.e., 

R d = 4A (ACDM) , (3.2) 

the contribution from Rd is negligible because A rj Hq w (10~ 33 eV^) 2 , where Hq is the 
present Hubble parameter. 



B. Exponential Gravity 



The exponential gravity has been studied intensively in the literature 22j, |47H54|. The 
form of f(R) in the exponential gravity model is given by 



f(R) = R - PRs (1 - e~ R / Rs ) 



(3.3) 



where R$ is the characteristic curvature scale and is a model parameter. The viable 
conditions are satisfied when (3 > 1 and R$ > 



47 



. The feature is that it is free from 



the fine tuning problem and it has only one parameter more than the ACDM model. 

Now we will first investigate the de-Sitter curvature R d in the model. The first and second 
derivatives of f(R) with respect to R are 



f{R) = l-Pe- R l Ra and f" (R) = -^-e~ R ' Rs . 
According to the condition for the de-Sitter curvature (12.41) . R d satisfies 

(1 - Pe- Rd/Rs ) R d = 2R d - 2(3R S (l - e - Rd/Rs ) . 



(3.4) 



(3.5) 



Defining x = R d /Rs, Eq. (13. 5p becomes 



x = 2(3- f3e~ x [x + 2) 



(3.6) 



The numerical solutions for this equation are shown in Fig. [TJ The factor e~ x (x + 2) 
decreases very fast when (3 > 1, which is generally required by the viable condition for the 
exponential gravity. Therefore, we can obtain the asymptotic solution of x for a large j3: 



x = 2(3 for > 1. 



(3.7) 



9 




Figure 1. Numerical solution of x = R^/Rs versus the model parameter /3 in the exponential 
gravity model, where the solid line indicates the exact solution of Eq. (|3.6p and the dashed line 
presents the approximated solution x = 2/3. 




Figure 2. m 2 s / R$ versus /3 in the region of /3 = to 2 (left panel) and /3 = 1 to 3.5 (right panel) in 
the exponential gravity model, where the solid lines indicate the numerical solution of Eq. (13. 8p . 
m 2 s /Rs is essentially zero when /3 = 1, and the dashed lines show the approximated solution x = 2/3. 



From Eq. (I2.16p . we derive the mass squared of the scalar mode in the exponential gravity 



as 



m 



The numerical solutions for m? s /Rs are presented in Fig. [2j Since in the large curvature 
regime R/Rs 3> 1, the theory will recover the cosmological constant model, R$ is roughly 
inverse proportional to (3 in the way that 



@R S = 2A = 9.94 x l(r 6( W 2 



(3-9) 
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Figure 3. m 2 versus {3 in the region of /3 = 1 to 2 (left panel) and /3 = 1 to 5 in log scale (right 
panel) in the exponential gravity model, where the solid lines present the approximated numerical 
solution of m 2 obtained by Eq. (|3.10p and the dots in the right panel show the exact value of m 2 
presented in Table [TJ 



Table I. Numerical results of the scalar mode mass m s in vacuum with respect to different /3 in the 
exponential gravity model. 



p 


h 


..ini 

Vh 


0° 


R s (10~ 6( W 2 ) 


m s (10~ 33 eF) 


4 


0.7050 


2.618 


0.2761 


2.452 


24.36 


3 


0.7059 


2.609 


0.2758 


3.263 


11.39 


2 


0.7103 


2.558 


0.2738 


4.824 


5.069 


1.27 


0.7194 


2.45 


0.2701 


7.39 


1.86 



with the value of A obtained from WMAP 7 (55fl, SDSS 7 56] and SCP Union2 observa- 



tions 



57] . Eq. (13. 8p then can be approximated as 

2A (\ 



3/3 \(3 6 



x 



1 



(3.10) 



In Fig. [31 we depict the result of the mass squared m 2 s versus (3 by using Eq. (13.1 OP . We also 
calculate the exact m s without any approximation, and the results of /3 = 1.27, 2, 3 and 4 
are shown in Tabled where we have used the values of Rs obtained from our previous result 
in Ref. J|33| under the constraints of WMAP 7, SDSS 7 and SCP Union 2 measurements. 
For ^> 1, the mass squared m 2 s becomes 

3 U 



(3.11) 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

X 



Figure 4. Numerical results of the de Sitter curvature Rd divided by R c versus A with n = 1 (solid), 
2 (dashed), 3 (dot-dashed) and 4 (long-dashing), respectively, in the Starobinsky model, where the 
dotted line presents the approximation solutions x = 2 A when A 3> 1. 

When (3 is O (1), m s is around 10~ 33 eV r . However, the cosmological observations do not 
give any significant upper bound on /3. Thus, m s could be arbitrary large in this case. As 
(3 — y oo, corresponding to the ACDM model with m s — > oo, the scalar mode of gravitational 
wave vanishes. 



C. Starobinsky Model 



In Ref. |36[, Starobinsky proposed the following f(R) form: 



f(R) = R-\R c [l-[l + ^) ), (3.12) 

where R c is roughly the present cosmological density and A and n are positive model parame- 
ters. From the solar system constraint and the bound on the violation of the equivalence prin- 



ciple, one gets n > 0.9 |2l|. In Fig. HJ we present the vacuum curvature Rd at the de Sitter 
stages obtained by Eq. (12. 4p . We can see that R^jRc — 2A when A ^> 1. Since for R S> R c , 
the model will restore the ACDM model, we have XR C ~ 2A and R d ~ 4A = 1.99 x 10~ 6f W 2 
when A> 1. 

In Fig. |5l we depict the numerical solutions of the mass squared of the scalar mode m 2 s 
derived from Eq. (I2.16p . 
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Figure 5. m 2 /R c versus A in vacuum in the Starobinsky model. 



IV. GRAVITATIONAL WAVE IN INNER GALAXY 

A scalar mode of gravitational wave with effective zero or non-zero mass in f(R) is a 
very different prediction from the ordinary GR. However, we will show that in the presence 
of matter density, the scalar mode might not be able to exist in these viable f(R) models. 
Consider a scalar mode of gravitational wave propagating within our Galaxy halo. The local 
homogeneous density of dark matter and baryonic matter is roughly p « 10~ 24 g/cm 3 . If we 
take this matter density into our analysis, it will give a large contribution to the background 
curvature compared to the vacuum de Sitter curvature. (The ratio of the matter density p 
to de Sitter curvature Rd is about n 2 p/Rd ~ K 2 p/AA m 10 5 .) In this case, the condition for 
the background curvature R f |2.4[) should be modified as 



where Rq is the background curvature with matter. Note that for viable f(R) models, the 
solutions to Eq. (14. II) can be approximated as Rq ~ n 2 p at the high curvature regime. 
In the case of the exponential gravity (13. 3ft . Eq. (14.1ft gives 



where x = Rq/Rs and r = n 2 p/Rs are the ratios of the background curvature and matter 
density to R$, respectively. Since 0Rs — 2A from ( 13. 9ft . we find that the solution of Eq. 
( 14. 2 p is extremely large, 



f'(Ro)R = 2f(R ) 




(4.1) 



x = 2/3 + r - (3e~ x (x + 2), 



(4.2) 




2 x 10 5 ^, 



(4.3) 



R d /2f3 
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Figure 6. m 2 versus n in the Starobinsky model with matter density p = 10 24 g/cm 3 and XR C = 2A. 

which just leads to Rq ~ k, 2 p. Thus, in the exponential gravity, the mass of the scalar mode 
will become an extreme in the galaxy region: 



The corresponding cutoff frequency uj m is also infinite. As a result, it is almost impossible to 
detect this scalar mode within our Galaxy under the exponential gravity scenario. Moreover, 
for any source that is massive enough to generate gravitational waves, we expect them to 
lay in the region with density higher than the baryonic/dark matter density 10~ 24 g/cm 3 . 
Therefore, the scalar mode will not have the chance to propagate from the source in the 
exponential gravity. 

In the case of the Starobinsky model (13. 12ft . the situation is quite different. The scalar 
mode of gravitational wave can have a light mass in the galaxy region. The minimum bound 
of the scalar mode mass is m s > 10 eV when p = 10~ 2A g/cm 3 . The corresponding cutoff 
frequency is quite small f m > 10 -9 Hz. This feature will allow the propagation of the scalar 
mode inside the galaxy. Hence, detecting the scalar mode in the Starobinsky model will be 
possible. In Fig. El we depict the mass squared versus the model parameter n with different 
fixed values of A, where we have used XR C = 2 A. The scalar mode can still be very heavy 
when the index n goes to a large value, but the mass dependence on the parameter A is not 
quite significant. 




(4.4) 
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V. CONCLUSIONS 



We have discussed gravitational waves in viable f(R) theories. Using the weak field ap- 
proximation on the field equation, we have confirmed that f(R) will give an extra massive 
scalar mode besides the ordinary tensor mode in the standard GR. We have explicitly in- 
vestigated the situations of the extra scalar mode of gravitational wave in the exponential 
gravity and Starobinsky models of the viable f(R) gravity theories. In vacuum, we have 
shown that the typical mass squared of the scalar mode is in the order of the de-Sitter 
curvature m 2 s ~ R^ « 10 _66 eV 2 in both models. 

However, in the galaxy region, the situations will be different if we consider the back- 
ground curvature by the presence of the galactic density. The small matter density p = 
10 -24 g I cm 3 is still 10 5 larger than the de-Sitter curvature Rj in both models. In the expo- 
nential gravity, since the mass of the scalar mode in galaxy is undetectable large, it will be 
unable to measure a gravitational wave. On the other hand, in the Starobinsky model, the 
mass can be much smaller with its lower bound in galaxy being about 10~ 24 eV (or 10~ 9 Hz). 
Therefore, it is possible to observe the scalar mode of gravitational wave in the Starobinsky 
scenario if there is an astrophysical source which generates this scalar mode. The extreme 
difference of the scalar mode masses also makes the distinction of f(R) theories from ACDM 
possible by probing this scalar mode of gravitational wave. 

Recently, there is an underway space-base d gr avitational wave probing experiment, the 
Laser Interferometer Space Antenna (LISA) 581 H which is a proposed joint mission of the 
European Space Agency (ESA) and NASA. It will measure the low-frequency band (10~ 5 to 
1 Hz) of gravitational waves with high signal-to-noise ratio. The gravitational sources within 



this band 



59] include supermassive black holes, intermediate-mass black holes, extreme- 



mass-ratio black hole inspirals, galactic compact binaries and some primordial gravitational 
wave sources. It is possible that these sources also generate the scalar mode of gravitational 
wave. Since LISA will be located far from the Earth and other gravitational sources, the 
background curvature of it is very low compared to the ground-based experiments, which 
allows the propagation of the scalar mode in some viable f{R). As a result, LISA has a great 
chance to direct detect not only the ordinary gravitational wave but also the clues of the 
deviation from Einstein's GR by analyzing the scalar mode behavior of gravitational wave. 



1 Laser Interferometer Space Antenna, http://sci.esa.int/lisa 
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Moreover, if the scalar mode becomes observable, many interesting features in the viable 



f(R) gravity mode 
as ASTROD-GW 



12 



60| will appear. We note that other gravitational wave probes, such 
6l| with the sensitivity in the 10~ 7 — 10 _1 Hz band, may also detect the 



extra scalar mode. It is clear that an observation of the scalar mode with a frequency larger 
than 10~ 8 Hz by LISA or ASTROD-GW would be associated with the Starobinsky model 
and is rules out the exponential one. Finally, we remark that the scalar mode in a viable 
f(R) cannot be observed by the ground gravitational searches due to the large background 
curvature. 
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